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Abstract .  This  paper  is  concerned  with  small  deformations  of  a  class  of 
laminated  layered  composite  plates  and  a  class  of  laminated  layered 
composite  rods  utilizing  a  direct  formulation  of  thermomechanical  theories 
of  shells  and  rods.  The  paper  is  arranged  in  two  parts,  namely  Part  A  (for 
laminated  composite  plates)  and  Part  B  (for  laminated  composite  rods) ,  and 
can  be  read  independently  of  each  other.  In  each  part,  after  providing 
the  main  ingredients  of  a  direct  formulation,  fairly  general  linear 
theories  are  developed  which  include  both  dynamical  and  thermal  effects. 
Some  applications  to  sandwich  plates  and  layered  rods  are  discussed  in 


detail . 


1. 


Introduction 


Among  the  many  types  of  composite  materials  which  have  been  studied, 
layered  media  have  received  much  attention.  One  group  of  such 
media  consists  of  a  large  number  of  alternating  plane  parallel  layers  of 
two  homogeneous  and  usually  isotropic,  elastic  materials.  Three  main 
approaches  have  been  used  to  construct  theories  to  describe  their  mechanical 
behavior.  One,  termed  the  effective  modulus  theory,  consists  in  replacing 
the  composite  by  a  homogeneous  anisotropic  medium  whose  material  constants 
are  determined  in  terms  of  the  geometry  and  material  properties  of  the 
constituents  of  the  composite.  Ihe  second,  called  the  effective  stiff¬ 
ness  theory,  yields  more  detailed  information  about  the  composite  than  the 
effective  modulus  theory;  and  is  expounded,  for  example,  by  Sun,  Achenbach 
and  Herrmann  (1968)  with  extensions  by  Drumheller  and  Bedford  (1973) .  This 
theory  also  has  some  connection  with  microstructure  theories  as  shown  by 
Herrmann  and  Achenbach  (1968)  and  Achenbach  (1975).  The  third,  discussed 
by  Stem  and  Bedford  (1972)  is  based  on  a  multi-continuum  theory  for 
closely  coupled  mixtures. 

Another  class  of  layered  media  consists  of  multilayered  shells  and 
plates.  The  earliest  theory  used  for  these  seems  to  have  been  of  the 
effective  modulus  type  for  elastic  composites,  which  gives  useful,  but 
limited  information  about  the  behavior  of  the  composites.  A  different 
approach,  used  by  Malcolm  and  Glockner  (1976)  for  a  three  layered  iso¬ 
thermal,  elastostatic  sandwich  shell  and  by  Epstein  and  Glockner  (1977) 
for  an  n  layered  elastostatic  shell  involves  a  single  surface  and  a 
single  director,  or  n  directors,  respectively,  one  associated  with  each 
layer  of  the  sandwich,  together  with  the  use  of  variational  procedures. 
Another  type  of  theory  given  by  Epstein  and  Glockner  (1979)  models  the 
elastic  shell  by  a  surface  and  n  directors  associated  with  the  surface. 
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and  again  employs  variational  methods.  Here  there  is  considerable  dif¬ 
ficulty  in  identifying  the  elastic  coefficients  of  the  sandwich  shell  or 
plate  in  order  to  render  this  theory  useful  for  specific  problems. 

The  effective  modulus  or  classical  theory  for  laminates  consisting  of 
linear  elastic  material,  has  been  extended  to  a  higher  order  theory  by  Lo, 
Christensen  and  Wu  (1977,1978)  and  has  been  discussed  also  by  Christensen 
(1979) .  This  theory  utilizes  one  expansion  for  displacements  throughout 
the  laminates,  whatever  the  number  of  laminates,  together  with  a  variational 
procedure,  and  has  been  derived  for  statical  problems  under  isothermal 
conditions.  Displacements  and  stresses  are  continuous  at  the  interfaces 
and  the  theory  yields  more  detailed  information  than  the  classical  theory. 

A  different  higher  order  theory  given  by  Pagano  (1978),  again  is  limited  to 
statical  problems  under  isothermal  conditions.  Pagano  (1978)  assumes  stress 
distributions  for  each  laminate  which  satisfy  appropriate  surface  conditions 
at  interfaces,  and  derives  the  displacement  fields  and  equations  of  equi¬ 
librium  with  the  help  of  a  variational  procedure. 

With  the  use  of  an  expansion  procedure,  Hegemier  and  Bache  (1973) 
have  discussed  special  types  of  one-dimensional  motions  in  a  plate  with 
isotropic  layers.  Similarly,  Fourier  integral  methods  were  used  earlier 
by  Buffer  (1971)  for  infinite  plates  with  isotropic  layers. 

In  the  present  paper  a  different  approach  from  those  noted  above  is 
used  to  obtain  a  linear  theory  of  elastic  sandwich  plates,  one  which  is 
available  for  dynamical  problems  and  which  allows  for  thermal  effects. 

The  method  used  can  also  be  applied  to  sandwich  type  shells  and  to  materials 
which  are  not  necessarily  elastic.  The  theory  of  plates  and  shells  for 
single  phase  media  has  a  long  history  —  see,  e.g.,  Naghdi  (1972).  In 
recent  years  considerable  attention  has  been  given  to  direct  two-dimensional 
theories  from  a  thermodynamical  point  of  view  where  temperature  effects 
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along  the  major  directions  of  the  plate  or  shell  are  considered.  The  full 
thermomechanical  theory  for  plates  or  shells  in  which  temperature  effects 
across  the  thickness,  as  well  as  along  the  major  surfaces,  are  incorporated, 
has  been  given  by  Green  and  Naghdi  (1979a).  In  particular,  the  linear 
thermoelastic  theory  for  a  plate  is  studied  in  detail  (Green  and  Naghdi 
1979a)  and  all  the  thermoelastic  coefficients  for  an  isotropic  plate  were 
found  in  terms  of  the  three-dimensional  coefficients.  This  theory,  extended 
here  to  orthotropic  plates,  is  immediately  available  for  each  individual 
member  of  a  sandwich  type  plate  (or  shell),  together  with  appropriate 
conditions  to  ensure  continuity  in  displacements  and  surface  forces  at 
interfaces,  or  appropriate  specified  values  for  surface  forces  or  displace¬ 
ments  if  other  conditions  prevail  at  interfaces.  The  final  form  of  the 
theory  has  some  similarities  to  that  given  by  Pagano  (1978)  when  the 
results  are  specialized  to  the  isothermal  static  case.  It  does  not, 
however,  yield  immediate  values  for  the  actual  stress  distribution  in 
each  laminate,  but  is  concerned  with  force  resultants  and  higher  order 
couple  resultants.  These,  if  desired,  can  be  used  to  obtain  approximate 
values  for  stresses  by  a  procedure  similar  to  that  of  Pagano  (1978). 

In  a  similar  way,  the  complete  direct  formulation  of  a  thermoelastic 
theory  of  a  rod  with  rectangular  cross-section,  which  has  been  developed 
by  Green  £t  al_.  (1974)  based  on  earlier  work  on  the  subject  (Green  and  Laws 
1966,1973),  is  employed  in  Part  B  (sections  10-14)  to  study  the  behavior 
of  laminated  composite  beams  in  which  each  layer  of  the  composite  is  a  plane 
rod  (or  a  beam)  of  rectangular  cross-section  with  different  orthotropic 
properties.  Although  our  main  developments  for  composite  rods  are  limited 
to  a  linear  isothermal  theory,  the  same  approach  together  with  the  use  of 
recent  results  for  nonisothermal  rods  (Green  and  Naghdi  1979b)  can  be  used 
to  account  for  thermal  effects  and  a  more  general  formulation  of  composites 
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in  the  presence  of  finite  deformation.  A  different  approach  from  that 
given  here  for  rods  has  been  used  by  Dokmeci  (1973)  who  considered  a 
nonlinear  isothermal  theory  of  elastic,  composite  beams  based  on  an 
expansion  for  displacements  in  each  layer  together  with  a  variational 
approach  for  obtaining  equations  of  motion.  Dokmeci  (1973)  obtains  his 
constitutive  equations  for  force  resultants  from  the  energy  equation  in 
the  three-dimensional  theory  but  this  portion  of  his  work  suggests  some 
doubt  about  the  values  of  some  of  his  constitutive  coefficients  especially 
those  for  flexure. 

Specifically,  the  contents  of  the  paper  are  arranged  in  two  parts  and 
can  be  read  independently  of  each  other:  One  part  (Part  A)  is  concerned 
with  the  linear  theory  of  laminated  layered  composite  plates  and  the  other 
(Part  B)  is  devoted  to  the  linear  theory  of  laminated  layered  composite 
rods.  In  Part  A  (sections  2-9),  first  a  concise  summary  of  the  basic 
equations  of  the  linear  theory  of  orthotropic  elastic  plates,  constructed 
by  a  direct  approach  via  the  linearized  theory  of  a  Cosserat  surface  with 
a  single  director,  is  given  and  this  is  followed  (section  3)  by  a  cor¬ 
responding  summary  of  a  restricted  linear  theory  of  plates  by  direct 
approach.  Next,  the  developments  in  sections  2  and  3  are  used  to  formulate 
a  linear  theory  for  a  laminated  plate  containing  N  laminates  each  consisting 
of  a  plate  of  homogeneous  orthotropic  material  of  constant  density  and 
thickness  and  each  at  constant  temperature  along  with  appropriate 
continuity  conditions  for  the  tractions  and  the  displacements  at  the 
interfaces. 

The  results  in  sections  4-6  are  specialized  to  composites  with  three 
laminates  and  are  applied  to  a  number  of  cases,  which  include  the  thermo¬ 
static  problems  of  the  torsion  of  a  three  layer  laminated  rectangular  plate 
in  the  context  of  thermo -statics,  the  study  of  harmonic  waves  and  the  case 
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of  periodic  force  resultant  distribution  in  a  stratified  medium  consisting 
of  a  large  number  of  alternating  parallel  layers  of  two  homogeneous 
orthotropic  materials. 

In  Part  B  (sections  10-14),  first  a  concise  summary  of  the  basic 
equations  of  the  linear  theory  of  orthotropic  elastic  rods,  constructed 
by  a  direct  approach  via  the  linearized  theory  of  a  Cosserat  rod,  is 
given.  These  are  then  applied  (sections  12-14)  to  some  simple  static  solutions 
for  composite  rods  and  harmonic  wave  propagation  along  the  composite  rod. 


Part  A.  Laminated  Layered  Composite  Plates 


2.  Basic  equations  of  the  linear  theory  of  orthotropic  elastic  plate 

★ 

Consider  a  plate  of  constant  thickness  h,  constant  density  p  and  at 
constant  temperature  9,  bounded  by  the  planes  z  =  ±^h,  where  x^  =  (x,y,z)  are 
rectangular  Cartesian  coordinates,  with  unit  vectors  e^  along  the  x^-axes. 

The  plate  has  homogeneous  orthotropic  symmetry  with  respect  to  the  z-axis 
and  two  axes  orthogonal  to  the  z-axis.  The  basic  equations  for  the  deforma¬ 
tion  of  such  a  plate  have  been  given  in  a  number  of  previous  papers  and,  in 
particular,  in  the  monograph  by  Naghdi  (1972)  which  also  contains  values  of 
the  mechanical  constitutive  coefficients  for  a  linear  isotropic,  elastic 
plate  (Naghdi  1972,  Ch.  E) .  The  full  thermodynamical  coefficients  for  an 
isotropic,  elastic  plate  have  been  given  by  Green  and  Naghdi  (1979a).  The 
main  equations  and  constitutive  coefficients  are  summarized  below  in  Cartesian 
tensor  notation  with  Greek  indices  taking  the  values  1,2.  Partial  dif¬ 
ferentiation  with  respect  to  xa  is  denoted  by  (  ) ,a  and  a  superposed  dot 
denotes  partial  differentiation  with  respect  to  time. 

We  recall  that  a  Cosserat  surface  C  comprises  a  material  surface  S  and 
a  director  assigned  to  every  point  of  the  material  surface.  In  the  reference 
configuration,  let  the  material  surface  be  a  plane  .SR  which  we  identify  with 
the  middle  plane  z = o  of  the  plate  described  in  the  preceding  paragraph. 

Let  u  and  6  denote,  respectively,  the  infinitesimal  displacement  vector  and 
the  infinitesimal  director  displacement  of  ,SR  (or  the  plane  z * o) .  Then, 
referred  to  the  basis  e^,  u  and  6  may  be  expressed  as 

u  =  u  e  +  u,e,  ,  6  =  5  e  +6,e,  (2.1) 

~  cua  3~3  ~  a~a  3~3 
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and  the  strain  measures  are 


ea3  =  3*(ua,6  +  uB,a)  ’  p3a  =  63,a 


o  -  =  <5  „  v  =  5  +  u, 

MaB  a, 3  Ta  a  3, a 


(2.2) 


We  denote  the  temperature  change  in  the  plane  SD  by  6  and  the  temperature 
change  along  the  thickness  of  the  plate  by  <p. 

The  contact  force  vector  N,  the  contact  director  force  vector  M,  the 
internal  force  director  m  and  the  entropy  fluxes  k,Z  are  given  by 


H  =  (NaB?3  +  Va?3)va  *  *?  =  (MaBS(3  +  Ma3!3)va  * 


m  =  V  e  +  V,e_  ,  k  =  p  v  ,  Z  =  p1v 
aua  3-3  a  ’  a 


where  v  =  va^a  is  the  outward  unit  normal  to  a  curve  in  the  (x,y)  plane  in  a 
specified  direction.  The  force  vectors  t,„t  on  the  major  surfaces  z  =  ±^h, 
respectively,  in  terms  of  their  components  are  given  by 


(2.3) 


t  =  t,e  +  t,e, 
u~  u  ct~a  u  3~3 


ll  =  Jlta!a  +  Jlt3?3 


(2.4) 


and  yk  ,^k  are  the  entropy  flux  vectors,  respectively,  at  these  surfaces. 
Then  the  field  equations  for  forces  and  entropy  fluxes,  in  which  body  forces 
and  entropy  sources  are  zero,  are 


NaB,a  +  ut6  +  ^3  =  pU3  ’ 


Ma3,a-V*  ^uVfcV  =  pa63  * 


(2.5) 


*  * 


-p  -  k  -  „k  =  pn  . 
*a,a  u  z 
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(2.6) 


Mo8.a-V*  h(uVtV  =  P“5B  * 


VaVsVs  =  pU3  * 

-Pa,a  +  pCl‘%  =  pT\  ’ 


where  n.Hj  are  entropy  densities  and  ^  is  an  internal  rate  of  production  of 
entropy. 

The  desired  thermoelastic  constitutive  equations  for  an  orthotropic 
plate  can  be  derived  from  a  Helmholtz  free  energy  function  ip  of  the  form 

2P*  =  <WXW*2%V3*A<53>2  *  EXUB3AP3u 

*  DxuDaepxu  *  Vx\ 

*  2B8,8  .  .  D62 

*  "xu^xu*  *  ^  • 

where  for  compactness  we  have  designated  the  usual  Cartesian  tensor 

/vQ  I’ufl 

constitutive  coefficients  A  ,D  as  A,  ,D,  defined  by 

aSAy  agAy  Ay  Ay  1 

A f  =  5a558nAr  .  ,  of  =  s'VV  , 

Ay  5nAy  ’  Ay  £nAy 

611  =  522  =  1  ,  512  =  621  =  0  . 


Also,  the  Cartesian  tensor  components  of  the  thermoelastic  coefficients 
satisfy  the  symmetries 


AagAy  “  AgaAy  "  AagyA  “  AAyag  *  AAy  AyA  * 


DagAy  =  DAyag  '  °Ay  DyA  ’  BAy  =  ByA  1 


The  corresponding  constitutive  equations  are 
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(2.13) 


10. 


33  33 

C33  “  oC33 


C33  =  oC33  ’  d33  =  od33  ' 


(2.18) 


In  order  to  specify  the  constitutive  coefficients  in  (2.10)  and  (2.11) 
in  terms  of  the  three-dimensional  coefficients  (2.16)  it  is  simplest  to  use 
axes  which  coincide  with  the  coordinate  axes  x,y  and  relate  simple  exact 
solutions  of  three-dimensional  elasticity  with  solutions  of  the  equations 
(2.S),  (2.6),  (2.10),  (2.11).  Coefficients  with  respect  to  general  axes 
in  the  x,y  plane  may  then  be  found  by  appropriate  tensor  transformations. 
The  results  are 


and 


Aa^  =  hca®  A  =  hc^  A  =  hc^ 

%  ncXy  *  aB  ncaB  ’  A  nc33  • 

Ea6=  6h  6aACB3  *  P  =  P  h  »  Bag  =  ’  hcag  *  B  =  "  hc33 


D  =  -  p  ch/9  ,  aag  *  hd^/0  » 


511  =  522  *  1  ’  S12  “  621  =  ° 


aB  33 

aB  _  J_  h3  aB  .  C33%, 

%  "  12  n  33  ;  ’ 

C33 

n  o',  r  A3  1  v2 

aB  =  Bh6aAcB3  ’  a  "  12  h  ’ 


(2.19) 


FaB  =  '  12  h  CcaB 


33 

aBC33 

33 

C33 


)  , 


(2.20) 


F  S_±h\els.+  !j3Cj3) 

h  12  n  1  5-  33  J  * 

9  C33 

,  h  .  .  1  h3  , 

b  =  g  d33  ’  baB  =  12  q  da6  ‘ 


An  approximate  comparison  with  the  torsion  problem  for  a  plate,  or  an  exact 
comparison  with  one  of  the  flexural  modes  of  wave  propagation,  leads  to  the 
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—  2  — 

choice  0  =  tt  / 12  for  3.  This  is  approximately  5/6,  the  value  specified  in 
previous  papers.  The  value  of  3  is  much  more  dependent  on  the  particular 
problem  chosen  as  a  basis  of  comparison.  Here  we  choose  3  =  as  an 
previous  papers. 

For  one  purpose  below  we  require  equations  of  the  form  (2.12)  for  a 
plate  which  has  symmetry  of  structure  at  any  point  only  with  respect  to 
directions  normal  to  the  middle  surface,  i.e.,  in  the  ^-direction,  but 
does  not  necessarily  have  symmetry  with  respect  to  the  x  and  y  directions. 
The  non- zero  coefficients  in  (2.12)  are  then  the  set  given  by  (2.16), 
(2.13)  and  (2.14). 
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Summary  of  a  restricted  linear  theory  of  plates 

In  addition  to  the  theory  summarized  in  section  2,  in  subsequent 

sections  we  wish  to  consider  also  the  use  of  a  system  of  equations  which 

correspond  to  those  of  a  classical  (Poisson- Kirchhoff)  linear  theory  for  a 

plate. having  symmetry  with  respect  to  directions  which  are  normal  to  the 

plate  but  which  is  not  necessarily  orthotropic.  These  equations  may  be 

obtained  from  the  present  equations  by  a  limiting  or  approximate  procedure. 

Alternatively,  they  may  be  derived  from  linearization  of  the  restricted 

theory  discussed  by  Naghdi  (1972,  Secs.  10,15),  but  extended  here  to  account 

for  full  thermal  effects  and  expressed  in  a  slightly  more  general  form.  The 

plane  faces  of  the  plate  are  given  by  z  =  a,  z  =  a  +  h  where  a  is  constant. 

Again,  we  identify  the  material  surface  <S  in  the  reference  configuration 

K 

with  the  plane  z  =  0,  which  is  any  suitable  reference  plane  in  the 
plate.  In  the  restricted  theory  the  director  is  coincident  with  the  outward 
unit  normal  to  the  plane  SD  so  that  <5=>B,  3  =  3  e  .  B  =  -u_  .  while  the 

infinitesimal  displacement  u  is  still  given  by  (2.1)^  The  relevant  field 
equations  in  the  restricted  linear  theory  are 


Na0  ,a  +  utB  +  «,tB  =  pUg  ’  u*3  '  J,t3  =  0  * 


*  ★ 


-p  -  k  -  „k  =  pn  , 
Mx.a  u  l 


MaB.a' V(a*h)uViMt|S  '  °  ’ 


V  +  t_  +  .t.  *  pu_  , 
a, a  u  3  l  3  3 


1  w  "  • 

-P  +  -  (a+h)  k  -  a  „k  =  pn,  , 

*a, a  KS1  ^  'u  l  y  1 


(3.1) 


where  for  a  plate  which  has  symmetry  only  with  respect  to  the  z-direction 
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and 


Nctg  =  Aa0eXy  +  A  agpXy  +  Bap9  +  Bi 


pn  =-  CVa6  +  Fipa6  +  D9  +  D/W  ' 


(3.2) 


P  =  -  a  ,  0  ,  -  a <f>  .. 
*a  aA  ,A  aA  ,A 


M  .  =  DApp.  +  A/APe.  +  F  +  F'  0  , 

a8  a0HAy  aS  Ay  a8  aB 


pril  =-  (Va8  +  Biea0  +  F0  +  D/0)  * 


Pa  *-baA^,A*baAe,A  ’  pCl  =  * b4>  ' 


(3.3) 


pa8  =  '  u3,a3  * 


Also  Va  is  not  now  determined  by  a  constitutive  equation.  The  various  coef¬ 
ficients  may  be  specified  in  terms  of  the  coefficients  (2.14)  by  studying 
simple  problems  which  have  exact  solutions  in  three  dimensions.  For  later 
convenience  we  assume  here  that  the  coefficients  (2.14)  are  functions  of  z. 
Thus 
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S*h  708  rhA, 

p  d*  •  *xu  I  CXnd  ’ 

3.  3 


7** 

Xp 


a+h 


C^zdz  ’ 


ra+h- 


oB 


Ca6dZ  ’  Bi  =  FaB  = 


D  = 


a+h 


C  dz  ,  D/  = 


a+h 
a 

ra+h 


C  Qzdz 
ci& 


C  zdz  , 


(•a+h 


“0  '  e 


Vz  ’  aaB  =  ba6=7a  daB 

a  o  J  a 


a+h 


Xp 


F  =  - 


a+h 

a 

a+h 


CXpz2(i2  *  FaB 


C  z2dz  ,  bft  =  ~" 
a  0 


a+h 

a 

a+h 


caez2dz 


d  .z  dz 
ci& 


-  i  ra+h 

b  =  —  d__dz  , 
9  U  " 


ciB  33 

raB  ciB  33CXp 
LXp  =  CXp  '  33 

C33 

33 

CaBC33 

CaB  =  CaB  *  33 

C33 


c  .  »*  s.  . 


9 


33 

'33 


zdz  , 


where 


4 .  Laminated  plate 

We  suppose  that  the  plate  contains  N  laminates  each  consisting  of  a  plate 

of  homogeneous  orthotropic  material  of  constant  density  and  thickness  and  each 

at  constant  temperature  0.  A  typical  laminate,  the  n1"*1,  has  constant  density 
★ 

p  and  thickness  h  and  lies  between  the  planes  z  =  H  ±  \  h  where  H  is 
n  n  r  n  n  n 

constant,  and  each  laminate  is  glued  to  its  neighbours.  Hence 


n+1 


-  H  =  %(h 


n+1 


hj  (n  =  1,2,...  ,N-1) 


(4.1) 


The  principal  directions  of  orthotropy  of  the  n^  laminate  are  at  angles  en 
to  the  x  and  y  axes;  without  loss  of  generality  one  value  of  en  (n  =  1,2, . . . ,N) 
can  be  taken  to  be  zero.  The  thermoelastic  coefficients  for  the  n^  laminate 
are  denoted  by  the  symbols  used  in  (2.12)  with  an  additional  prefix  n;  thus 
ncTj,  nc..,  nd^ ,  nc  and  the  relevant  coefficients  are  specified  by  (2.13), 
(2.14)  and  (2.16).  These  coefficients  may  be  expressed  in  terms  of  the 
principal  values  (2.15)  and  the  angle  by  the  usual  tensor  transformations. 

Each  laminate  is  assumed  to  be  governed  by  plate  equations  of  the  type 
given  in  section  2  with  constitutive  coefficients  specified  by  nA^, .  . .  ,nF, . . . , 
nb  q .  These,  in  turn,  are  given  in  terms  of  the  coefficients  in  (2.16)  by 
formulae  of  the  type  (2.19)  and  (2.20),  each  coefficient  having  the  prefix  n 
as  in  nc^j.  The  displacement  of  the  middle  plane  z  =  of  the  n^  laminate, 
together  with  its  associated  director  displacement,  are  denoted  by 


n  n  n 
"n  =  uci~a  u3~3 


6  = 
~n 


-n 
o  e 
a~a 


,n 

+  63!3 


(4.2) 


The  two  temperatures  associated  with  the  laminate  are  denoted  by  ©n , 4>n -  Since 
displacement  and  temperature  are  continuous  at  each  glued  face  of  the  laminate, 
it  follows  that 
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u  ,  -  k  h  .<5  =  u  %  h  <5 

~n+l  n+l~n+l  -n  n~n 


6  ,  -  %  h  ,<f>  ,  =  0  +  ^  h  <{> 

n+1  n+1  n+1  n  n  n 


,  (n-l,2,...,N-l)  . 


The  forces  acting  on  the  faces  z  1  H  ±  %  h  of  the  nt  laminate  are 

n  n 


n  .n  .n  .n'  ^n  ^n 

u~  =utaSa  +  ut3!3  ’  ll  =  lZc£a  *  Zt3~3 


(4.4) 


★n 

and  the  corresponding  entropy  fluxes  are  ^k  ,^k  .  Because  forces  at  the 
boundary  of  two  laminates  are  equal  and  opposite  and  entropy  fluxes  are 
continuous,  we  have 

n  n+1  „  .  *n  ,  *n+l  _  ..... 

t+.t  =0,  k  +.k  =0  (n  =  1,2, . . .  ,N-1) 

u~  JU  u  l 

Forces,  entropy  fluxes  and  temperatures  at  the  outer  boundaries  z  =  -  k  hj, 

z  =  H^  +  %  hN  of  the  plate  are  respectively 


li1’  *k#1,  ei-%h*i  and  u!N*  uk*N'  0N^  Vn  • 


Although  we  have  considered  only  surface  conditions  for  laminates 
attached  continuously  to  its  neighbors,  the  above  conditions  may  easily  be 
modified  to  allow  for  other  interfacial  conditions.  Force  and  director 
force  vectors  and  tensors,  entropy  fluxes,  vectors  and  densities,  etc.,  for 
the  nth  laminate  are  denoted  by  the  symbols  used  in  section  2  with  an  added 
index  n.  Thus,  e.g.,  we  have  , . . .  ,qn ,n" , . . .  ,p*n.  The  equations  of 

motion  and  entropy  balances  for  the  n^  laminate  are  of  the  same  form  as 
(2.5)  and  (2.6)  with  the  added  index  n,  so  that 
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VIn  *n  „  -n 

N  „  +  t„  +  *  P  U„ 

ctB.a  u  B  A  B  n  B 


m".  (t"-.t")  --paS"  , 

a3,a  3  nni  3  A3'  Kn  n  3 


(4.7) 


n  . *n  .  *n  „  *n 
-p  -  k  -  „k  =  p  n 
^a,a  u  A  n 


Vu'S-t1?’  -Vn'e  ■ 


,ji  n  n  -n 

v  +  t_  + „ t_  =  p  u_  , 

a, a  u  3  A3  Hn  3 


(4.8) 


In  _n  ,  ,  ,  *n  .  *n.  "n 

Pa,a  +  Pn?l"*  hn(uk  ‘  Ak  >  =  Pnnl 


for  n=l,2,...,N.  The  constitutive  equations  are 


,,n  n.Xu  n  n.  ..n  n.  „ 

“ce  ‘  Aae%  *  Ws  *  bm6„  ■ 


n  _  n  n  n»  ^ri  xi-»  « 

7  A  0  +  A  0,  *  B  9 
3  a0  ag  3  n 


Hi. 

Ma3  “  EaXP3X  * 


(4.9) 


n  rTin  n n  n  ^  \  n  n  a 

pnq  =-(B<S3+  Ba6ea6+  D  0J  ,  Pa  -  -  a^B^  . 


^  n  Xy  n  n„  ,  .xi  _  n  n 

aB  '  DaBpuX  Vn  *  a  "  DaXYX  ' 


pn  n*  .  n  /H_  n  n«  i  > 

pn?i  =  -  b  4>n  ,  PnnL  -  -  (  Fa6PaB  ♦  F  V  > 


(4.10) 


1°  ni.  A 

P«  =- 


aByn,X 


Equations  (4.3)  give  N-l  relations  between  the  2N  variables  u^.d^ 

(n  =  1 , . . . ,N)  so  that  some  restriction  is  placed  on  the  freedom  of  choice  of 
boundary  conditions  at  the  edges  of  the  plate.  To  determine  possible  choices 
of  boundary  conditions  we  consider  the  rate  of  work  of  external  forces  Nn 
and  external  director  forces  M11  which  is 


N 

R  *  Z  (Nn  •  un  +  NT  •  6r)  . 

n=l 


(4.11) 
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Now,  from  (4  2)  we  have 


U2  =  ^{^2^2  +  ^1^1^  +  ui  » 
n-1 

~n  =  ~1 +  ^(hn~n  + hl~l}  +  1  hr~r  0*»3f...,N)  , 


so  that  (4.11)  may  be  written  in  the  form 

N  .  N 

R  *  Uj  »  ZNn  +  61'(r*!ih1  E  N11) 
n=l~  ~  1  n=2 

+  5n  •  Qp  *  k  hNf) 

N-1  .  N 

tl  5  ■(Mnt4hNnih  E  Nr) 
n=2  ~n  n~  n  r=n+l~ 


(4.12) 


(4.13) 


The  expression  for  R  indicates  that  boundary  values  may  be  chosen  for 
or  their  coefficients  in  (4.13). 

Similarly,  in  view  of  (4.3),  some  restriction  is  placed  on  the  freedom 
of  choice  temperature  or  entropy  flux  boundary  conditions  at  the  edges  of 
the  plate.  The  total  heat  flux  at  an  edge 
N 

h  =  E  (kn0  +  ) 

n-1  n  n 

N  .  N  MM 

-  0X  E  kn  +  <fr  (SI +  k  hj  E  kn)  ♦  $  (JlN  hnkN) 
n=l  n=2  n 

N-1  N 

+  E  <J>  (An  +  H  h  kn  +  h  E  kr)  .  (4.14) 

n=2  n  n  n  r_n+^ 

Boundary  values  may  be  chosen  for  8^,<^, . . .  ,4>N  or  their  coefficients  in 
(4.14).  Alternatively,  N+l  boundary  conditions  of  a  suitable  type  may  be 
prescribed  at  each  boundary. 

The  above  theory  of  laminates  provides  values  for  >V^,  ^3  but  not 


As  is  customary,  through  appropriate  definitions  of  stress  resultants 
(see  Naghdi  1972,  Secs.  11-12),  the  quantities  ma^  be 

interpreted  as  membrane  force  resultants,  bending  moments  and  trans¬ 
verse  shear  resultants,  respectively. 
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for  stresses.  [In  fact,  any  plate  (or  shell)  theory,  whether  constructed  by 
a  direct  approach  or  derived  from  integration  of  the  three-dimensional 
equations  across  the  thickness  of  plate  (or  shell),  can  only  involve  forces 
or  resultants  and  not  the  components  of  the  stress  tensor.]  Information 
about  the  latter,  if  desired,  may  be  obtained  by  introduction  of  an 
additional  assumption.  In  particular,  in  the  context  of  the  present  paper, 
some  estimation  of  the  three-dimensional  stress  distribution  across  each 
laminate  may  be  obtained  by  a  procedure  similar  to  that  used  by  Pagano 
(1978).  Approximate  expressions  for  the  in-plane  stress  components  t^ 
for  the  n^  laminate  which  are  consistent  with  the  definitions  of 
as  stress  resultants  are 


.  n  1  ,.n 

tag  h  a$  +  . 
n  h 


12z  „n 
3  a@ 


th 


(4 


With  the  help  of  the  three-dimensional  equations  of  motion  for  the  n 
laminate  together  with  the  surface  conditions  (4.4)  and  equations  (4.7)  and 


(4.8),  we  may  obtain  the  following  expression  for  t^3  and  t^j,  namely 
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.  n  j,n  *  f  .  n  .  n»  z  *  n  .  n> 

g3  =  t3g  =  *^6  '  «,V  +  h”  CutB  +  S,V 

n 


+  ±(iz2-hh2n)[-\%  +  Hhj„t*-gt*))  , 


nSiB  l  SJ 


(4.15) 


n  n  1  f .  3  1  «  2  x  .  n 

t33  =  2  tg3,B  -  ^  (4Z 
1  n 


tn  1 

WB,8  "  O,0J 


.  82  ,2  ,  ,2,,,n  1  ,,  2  ,  ,2.,  n  , 

x  -  -4  h  )  vQ  „  -  (3z  -  4  h  )  (  t.  .  +  t.  _) 


V*B,B  2h 

»i  it 

n 

n 


n'  VB,ro,&J 


1  fry  2  «  ,  2.  f  n  . H\ 

.-J  (3!  \Huts-^1 

n 


z  ,  .n  ^.n. 

T-tuViV  • 

n 


(4.17) 


We  note  from  (4.16)  and  (4.17)  that 


fa h 


“  *6J**  *  1 


kh 


fah 


Ct33‘  Ztg3,B)dz  = 


(4.18) 


in  comformity  with  the  values  of  but  (4.16)  does  not  fit  in  with 

the  value  for  This  is  partly  due  to  the  nature  of  approximation  in 

the  special  assumption  (4.15). 
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5.  Effective  modulus  theory 


We  summarize  here  an  effective  modulus  theory  similar  to  the  classical 
theory  for  sandwich  plates  of  the  type  described  at  the  beginning  of  section  4. 
The  theory  is  dynamical  and  includes  thermal  effects.  The  composite  plate 
is  regarded  as  a  single  plate  with  symmetry  with  respect  to  directions  normal 
to  the  plate,  and  its  motion  is  governed  by  equations  (3.1)  to  (3.3).  The 
coefficients  in  these  equations  are  specified  in  terms  of  the  constant  coef¬ 
ficients  of  the  separate  laminates  by  using  the  formulae  (3.4). 

With  the  help  of  the  values  (2.19)  and  (2.20)  for  each  laminate,  the 
coefficients  in  (3.1)  to  (3.3)  have  the  values  listed  below: 


6. 


Simple  static  problems  for  a  laminated  plate 


We  consider  first  the  simple  homogeneous  extension  of  the  plate,  using 
the  theory  of  section  4.  Recalling  (4.2)  and  the  conditions  (4.3),  we  sup¬ 
pose  that  the  components  of  the  displacement  and  the  director  displacement 
vectors  and  temperatures  of  each  laminate  are  given  by 


Uj  =  Lx  +  My  ,  u”  =  Mx  +  Ny  ,  6^  =  0  ,  0^  =  0  ,  $n  =  0  , 

. N-‘>  ■ 

where  L,M,N, 0,^,6^  are  constants.  Then, 

ell  =  ®11  =  L  ’  e12  =  ®12  =  M  ’  e22  =  e22  =  N  ’ 

pnc  =  0  ,  p“  =0  ,  yn  =  0  . 
a3  3a  a 

With  the  help  of  (4.9)  and  (4.10),  it  can  be  seen  that  equations  (4.S), 
(4.7)  and  (4.8)  are  satisfied  provided 


t"  =  0  ,  0t"  =  0  ,  t™  =  0 

u  g  *  l  g  u  3 


*<3=0  *  uk*n  =  0  *  *kn=° 


for  n  =  1 ,2 , . . . ,N  and 


Vn  r*  n.  xn  n  « 

3  =  .  ,  A  53  =-  A^-  B  0n  . 


By  (6.4)2,  the  kinematical  quantities  6^  are  determined  in  terms  of  the 
constant  strain  components  =  eag  and  constant  temperature  0^  =  0. 
Then,  u^  may  be  an  arbitrary  constant  and  u^  is  given  by  (6.1)  for 
n*2 . N.  From  (6.4),  (4.9)  and  (4.10)  we  have 


(6.1) 


(6.2) 


(6.3) 


(6.4) 


n .  n .  n  * 

„n  ,n,\y  ag  a  y,,  ,n„  <2(3 

"as  "  <  Ac8  •  '.  '  ~)e>»  *  <  B«S  '  -^-,6  ■ 


nA  „nB  nBn. 

n  <-no  ag  ..  /H-  B  B,.  n  n 

pnn  =  -  (  Bag  '  n~  eag  "  (  D  ’  ~ >0  *  Pa  =  0  * 


(6 


Making  use  of  (5.1),  we  calculate  the  expressions  for  the  forces  N^, 
director  forces  Mag  with  respect  to  the  reference  surface  <SR  (or  z  =  0)  in 
section  4,  and  entropy  densities  for  the  whole  plate,  in  the  form 


N  _. 

S  .  =  n".  =  AApe.  +  B  0e 
ag  ,  ag  ag  Ap  ag 
n=l 


N 

pr,  =  I  pnnn  =  -  BageaB  -  D0  , 
n=  1 


n=i 


(6 


pnlsI,HnV  Bageag 
n=  1 


D;0 


Comparison  of  (6.6)  with  the  constitutive  relations  (3.2)  and  (3.3)  shows 
that  for  homogeneous  extension  of  the  composite  plate  the  overall  moduli  are 
the  same  as  those  used  in  the  effective  modulus  theory  of  section  5. 

Next,  we  consider  flexure  of  the  plate  with  the  reference  surface 
iS  (z =  0)  unstretched  and  the  entire  plate  remaining  unchanged  in 
temperature.  Again  recalling  the  continuity  condition  (4.2)  and  the 
relations  (4.1)  we  take  displacement  components  and  tenperatures  in  the 
forms 
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0 


(6 


<5^  =  Px  +  Qy  ,  =  Qx  *  Ry  »  ua  =  Hn5a 

Uj  =  -  %  Px2  -  Qxy  -  %  Ry“  +  Kn  ,  <J>n  =  0  ,  9^=0  , 

where  P,Q,R,  Kn,<|>  are  constants.  Also  fi”  are  constants  subject  to  the 
conditions 


k"*1  -  K"  .  H  .  H  (n  =  1.2 . N-l) 


(6 


Then, 


pll  =  pll  =  P 


P21  =  P12 

=  Q  » 

n 

n 

:ag  =  Hnpag 

>  Yi 

522  “  P22  =  R 


4"> 


(6 


Equations  (4.7)  and  (4.8)  are  satisfies  provided  (6.3)  and  (6.4)  hold,  where 
(6. 4) ^  determine  6^.  The  constants  Kn  can  then  be  found  from  (6.8)  with  K* 
being  chosen  arbitrarily.  From  (6.9),  (6.4),  (4.9)  and  (4.10)  we  have 


nA  nA 

.ji  ,,  ,n.Ay  a(3  Ay.  n  „  ,nD 

N  „  =  H  (  A  o - - - -)p,  ,  pn  =-H(B0 

ag  nv  ag  n  JKAy  ’  Kn  nv  ag 


nA  nB 

ag  V 

n .  pag  ’ 


,n  n_.Ay  n  n„ 

Ml.„  =  D^op,..  ,  p  n,  =-  F~p, 


ag  agHAy  ’  Kn'l 
From  (6.10)  and  (5.1)  we  then  obtain 


agMag  * 


(6 
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C6.ll) 


» 


where  and  are  computed  with  respect  to  the  reference  surface  >SR 
(or  z = 0)  in  section  4.  Comparison  of  (6.11)  with  the  constitutive 
relations  (3.2)  and  (3.3)  for  pure  homogeneous  bendir.j  in  which  eag  =  0> 

<(>  =  0,  0  =  0  shows  that  the  overall  bending  moduli  are  the  same  as  those 
used  in  the  effective  modulus  theory  of  section  5. 

In  the  foregoing  examples  the  overall  moduli  for  simple  extension  and 
flexure  of  the  composite  plate  are  the  same  as  those  found  from  the  effective 
modulus  theory,  but  this  is  not  generally  true  for  more  complex  problems. 

For  example,  if  the  plate  is  subjected  to  a  constant  temperature  difference 
between  its  major  surfaces,  then  the  overall  moduli  differ  from  those  for 
the  effective  modulus  theory.  Moreover,  for  more  complicated  problems  it 
is  difficult  to  complete  the  analysis  based  on  the  theory  of  section  4 
without  specifying  a  particular  value  for  the  number  N  of  laminates.  On 
the  other  hand,  it  is  no  more  difficult  to  use  the  effective  modulus  theory 
of  section  5  for  N  laminates  than  for  1  or  2  laminates.  For  many  purposes 
this  theory  may  be  adequate  but  sometimes  a  more  detailed  theory,  such  as 
that  in  section  4,  is  required.  For  example,  studies  of  simple  harmonic 
wave  propagation  by  the  effective  modulus  theory  fail  to  indicate  any 
dispersion  effects. 

The  theory  of  section  4  can  be  used  as  a  basis  for  a  numerical  discussion 
of  thermo-mechanical  response  in  a  composite  once  N  and  the  composition  of 
the  laminates  are  given.  In  the  next  two  sections,  we  limit  our  attention 
to  further  analysis  for  one  particular  composite  consisting  of  3  laminates 
of  orthotropic  material  whose  directions  of  orthotropy  coincide  with  the 
x,y  axes,  the  two  outer  laminates  having  the  same  thermo-elastic  moduli. 
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7. 


Composite  with  three  laminates 

In  the  theory  discussed  in  section  4,  we  now  restrict  the  axes  of  ortho- 
tropy  of  each  laminate  to  coincide  with  the  x,y  axes  so  that  non-zero  constitu¬ 
tive  coefficients  for  each  laminate  are  given  by  (2.15),  with  an  added 
prefix  n,  and  by  (2.13)  and  (2.14).  The  composite  consists  of  three 
laminates  in  which  the  outer  ones  have  the  same  mass  density,  thickness 
and  elastic  moduli.  We  choose  the  reference  plane  SD  (or  z » 0)  to  be  mid- 
way  between  the  faces  of  the  composite  plate.  Then, 


h3  =  hx  ,  H2  =  0  ,  H3  =  -Hj  =  %(h2+h2)  , 


★ 


★ 


Ap 


=  !AaB 
Xp 


etc. 


(7.1) 


Evidently,  for  3- layer  composite,  the  equation  of  motion  (4. 7)- (4.8)  and 
the  constitutive  equations  (4. 9) -(4. 10)  separate  into  two  groups:  One 
group  represents  the  extensional  motion  (or  the  stretching)  and  the  other 
characterizes  the  flexural  motion  (or  the  bending) .  In  order  to  record  the 
relevant  equations  in  an  economical  manner,  we  introduce  the  quantities 
ui’  Ui’  etc- »  defined  as  follows: 
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‘Pa 

»  Pa 
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Then,  with  reference  to  the  3-layer  composite,  the  first  group  of  equations 
appropriate  for  extensional  motions  are  given  by 


and  by 
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VaVsVsVlVs-0!^  • 
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(7.5) 
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The  second  group  of  equations  appropriate  for  flexural  motions  of  the 
3- layer  composite  plate  is 
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ua  =-h25a'^  hl6a 


U3  =  2u3  ■  hl53 


0  = 


-h2$2  -H 


(7.12) 


In  the  remainder  of  this  section  we  consider  an  example  in  the  context 
of  thermostatics,  namely  the  torsion  of  a  three  layer  laminated  rectangular 
plate.  Let  the  plate  be  bounded  by  the  planes  x=±^a,  y  =  ±  %  £,, 
z  =  ±(hj+%h2)  and  let  the  plate  be  subjected  to  the  action  of  couples 
about  the  y  axis  over  the  ends  y  -±H  l.  For  this  problem,  4^ ,<|)2 ,<j>3. 
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©1 »©2 *®3  are  zero  anc*  we  choose  the  following  forms  for  the  components 
of  displacements  and  director  displacements  which  are  compatible  with  (7.7) 
and  (7.12): 


1  u  3  u  2 

U1  =  HlTy  ’  ui  =  '  HiTy  »  u\ 


.2  x2'  „ 

fij  -  51  =  Ty  , 


Uj  *  -  h  T[h2A(x)  +hxB(x)]  ,  \  x[h2A(x)  +  h1B(x)]  ,  *  0 

<$2  =  ^2  =  tBM  *  Sl  =  tAW  > 


1  2  3 

u3  =  U3  =  U3  =  -  Txy 


12  3 

53  ■  53  ■  S3 


(7.13) 


where  t  is  the  twist  per  unit  length  and  A,B  are  functions  of  x  only.  It 

follows  that  u+,u,,6  ,6*  are  all  zero  and  relevant  equations  are  (7.8)  to 
Ot  o  Ct  j 

(7.12)  for  flexural  motions.  The  components  of  displacements  and  measures 
of  deformations  are 


u^  =  2H1xy  ,  u~2  =  -  x[h2A(x)  +hxB(x)]  ,  u*  =  ~  2xxy  , 

6^  =  2xy  ,  <5*  =  2xB(x)  ,  5”  =  0  , 

P12  +  P21  =  T('^x+i)  *  y2  =  t(:a"x)  »  y2  *  2tCB-x)  , 


*  i  rt_  ,*cid  «  n  ,  rdA  ni  n  ^ud  ,  ^ 

e12  ^VdT  1)+h2(dT  ^  *  P12  +  P21  =  2T(d?+1)  * 

the  remaining  components  of  strain  being  zero.  It  follows  from  (7.4), 
(7.6),  (7.9)  and  (7.11)  that  the  components  of  various  forces  and  director 
forces  are  given  by 


(7.14) 
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N12  *  N21  =  '  A12T^hltdx*  +h2Cdx*1^  ' 


(7.15) 


MU  '  H21  '  2  11  '  V2  '  2  V»-*>  ' 


3  1 

the  remaining  components  being  zero.  Given  that  the  forces  t.  ,  0t.  on  the 

u  1  X»  X 

faces  of  the  plate  are  zero,  the  equations  of  equilibrium  (7.3),  (7.5), 
(7.8)  and  (7.10)  are  satisfied  provided 

7  21  ^^12  12 

-#  -  i  w!  *  u*2>  ■ 0  •  -sr^Vi-0  ■ 


11  12  +  12  12 
—7——  -  V,  +  %  h. (  t, ♦  t,)  =  0  ,  t,  -  ,  tt  , 

dx  2  1  u  2  u  2  u  2  u  2 


(7.16) 


„1  *2  «-2  n 

t,  =  t.  =  t,  =  t_  =  0 

ul  ul  u3  u3 


Substituting  (7.15)  into  (7.16)  and  eliminating  t,  +  we  obtain 
the  following  differential  equations  for  A  and  B: 

b  +  d^B  _  d(A_x)  .  o  > 
d  xT  dx“ 

c  ~  +  -  e(B-x)  =  0  , 

dxz  dxz 


(7.17) 


where  the  coefficients  b,c,d,e  are  given  by 
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P12  'f  ^  h2  A12 

^  hlh2  A12 


1 _ 1 2  ,  v2  1.12 
P12  *  ^  hl  A12 

1  v  u  1.12 

*  hlh2  A12 


(7.18) 


d  =  .  .  1  12  *  6 
hlh2  A12 


,  u  1.12 

hlh2  A12 


We  need  to  supplement  the  differential  equations  (7.17)  with  suitable  boundary 
conditions.  To  this  end,  we  first  note  that  in  the  torsion  problem  under 
discussion  there  are  no  applied  forces  at  the  edges  x=±a.  Then,  in  view  of 
(4.13)  for  the  special  case  N=3  and  recalling  the  nonvanishing  expressions 
in  (7.15),  the  appropriate  edge  conditions  are 


M12  "  %  hlN12  =  ° 


,  -  h  h2N12  =  °  0n  x  =  1  a  • 


(7.19) 


The  above  conditions  are  satisfied  if 


~7“  +1  —  0  ,  ^-+1  =  0  on  x=±a 

dx  dx 


(7.20) 


From  (7.17)  and  (7.20),  it  follows  that 


A  =  x  +  L  sinh  nx  +  M  sinh  mx  , 


2  2 
_  bn  -d  ,  .  ,  bm  -d  ,,  .  , 

B  =  x  -  — j L  sinh  nx  -  — j —  M  sinh  mx  , 

n  m 


(7.21) 


where 


LdCm^-n2) cosh  na  +  2n [  (b+l)m2  -  d]  =  0  , 


Md(m2-n2) cosh  ma  -  2m[ (b+l)n2  -  d]  =  0  , 


(7.22) 


and  n,m  are  the  positive  roots  of  the  equation 


(bc-l)X4  -  (dc+be)X2  +  de  =  0 


(7.23) 


The  couple  at  the  edge  of  the  plate,  measured  per  unit  length,  is 
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given  by 


2  2 
(xej^  +  ye2)  *  N  +  e3  x  M 


+  (xe^  +  ye2  +  H^)  xfl'n^M1 


3  3 

+  (xe^  +  ye2  +  H^e^)  *  N  +  x  M 


Over  any  edge  y=  constant,  with  the  help  of  (7.15),  (7.24)  becomes 

SlXtxf*  V2)*52[M21*M21*H1N21'X(V2*V2)1  • 

This  yields  a  resultant  couple  tG  over  v = constant  about  the  y-axis  given  by 


tG  = 


{M21*M21*HlN21-*CV2*V2>}d*  • 

-a 


It  can  be  seen  from  (7.16)  that 


V2  *  V2  *  £  <M21  *  M21  *  HlN2l’  • 


so  that  with  the  help  of  this  result  and  the  edge  conditions  (7.19),  (7.26) 
becomes 


tG  =  2 


-a 


{M21  +M21  +HiN2l^dx 


With  the  help  of  (7.15),  (7.18),  (7.21)  and  (7.22),  it  follows  from  (7.28) 
that  the  torsional  rigidity  of  the  laminated  plate  is 


G  =  4h.  h_  1A;?[(b+c.2)a  ♦  {d  -  lb+l)m  H,Cl:bc)r l_^c1l]}tanh_ na 
1  1  u  nd(m  -n  ) 

(d  -  (b+l)n2){(l-bc)m2 + d(c+l)}tanh  ma  ■, 

2  2  J  • 
md(m  -n  j 

It  can  be  shown  that  when  h^-*-0  (7.29)  reduces  to  the  known  value 


i 


(7.24) 


(7.25) 


(7.26) 


(7.27) 


(7.28) 


(7.29) 
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(7.30) 
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3  U  ma  ’ 
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h,  (2  12  S  hl 
2  c, 


2  12  2  23 

for  a  plate  of  thickness  h2  and  moduli  c^2>  c23  computed  from  Cosserat 

plate  theory.  Similarly,  when  h2"*0,  (7.29)  reduces  to  the  expected  value 


,,  .3  1  12 

16a  h.  c. _  .  , 

1  12  ,,  tanh  ma, 

3  U  "  ma  1 


1  23 

_  -  JL. 

m  ‘  2h.  (1  12-1 
1  c12 


(7.31) 


1  12  1  23 

for  a  plate  of  thickness  2h1  and  moduli  c^ »  c2 3' 

When  the  width  2a  of  the  plate  is  large  compared  with  both  thicknesses 
h^.h^  of  the  laminates,  then  na  and  ma  will  be  large  and  the  torsional 
rigidity  (7.29)  reduces  approximately  to 


G  =  4h1h2a(b+c+2)1A^ 


(7.32) 


Also,  with  the  help  of  (7.18)  and  (5.1),  the  result  (7.32)  becomes 


—12 

G  =  8a  Dj2 


(7.33) 


This  is  the  same  value  which  would  be  obtained  by  using  the  effective  modulus 
theory  of  section  5.  However,  if  2a  is  not  large  compared  with  both  h^  or 
h0  the  effective  modulus  theory  is  not  sufficient  for  the  discussion  of 
torsion.  This  is  in  contrast  to  the  situation  for  pure  tension  and  flexure 
of  the  plate,  where  the  effective  modulus  theory  gives  satisfactory  values 


for  overall  moduli. 
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8.  Harmonic  waves  in  3-layered  plate 


We  consider  here  wave  propagation  in  a  3-layered  composite  plate  of 
the  type  discussed  in  section  7  utilizing  the  two  groups  of  equations  (7.3) 
to  (7.7)  and  (7.8)  to  (7.12).  We  limit  our  attention  to  the  propagation 
of  harmonic  waves,  under  isothermal  conditions,  along  the  x-direction  of  a 
composite  plate  whose  major  faces  are  free  of  forces  and  couples.  We 
assume  that  all  variables  contain  a  factor  exp  i(£x-cot),  where  are 
constants,  so  that  this  factor  can  be  removed  throughout  the  calculations 
that  follow.  The  group  of  equations  (7.3)  to  (7.7)  then  separate  into  two 
further  groups,  the  first  being  concerned  with  symmetric  SH  waves  and  the 
second  with  symmetric  plane  strain  waves.  For  SH  waves  we  have  non-zero 
values  for  u 2>®2'u2 


2  2  n  5  +  1  12  + 

N12  =  ^  A12^U2  ’  M12  =  ^  A12')u2  ’ 

m"  =  ie(1oJ^3<s;  .  v:  *  V-a;  ,  u!  =  2u^-%h .<$:  , 


2  2  1  2  2 
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12  u  2  u  2  2  2 


(8.1) 


+  A  -  A  -  -  pl“2u2  ' 


i5M‘l2  -  v2  +  *  hlAz  -  • 


A  non-zero  solution  of  equations  (8.1)  leads  to  the  dispersion  relation 


4  2  2  2  4 

(b-coa)e  -*»♦«»-  0  , 


(8.2) 
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where 
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e  -  %  D*hf C2Pjh2  -  p'hj)  , 

and  where  use  has  been  made  of  the  expressions  of  the  form  (2.19)-(2.20) . 

2 

The  dispersion  relation  (8.2)  yields  two  expressions  for  £  as  a  function 

2 

of  ai  ,  corresponding  to  two  different  modes  of  propagation  of  SH  waves.  In 
one  of  these 


2  d  2 
^  -*■  b  “  = 


p->h2  +  2p1h1 

~  - 12  ^ ,,  m 

h2  C12  2hl  C12 


(8.4) 


as  ai  -*0.  This  agrees  with  the  values  obtained  for  SH  waves  from  the 
effective  modulus  theory  of  section  5.  [Recall  that  there  is  no  dispersion 
relation  in  the  effective  modulus  theory.]  In  the  case  of  isotropic  layers 
the  limiting  value  (8.4)  as  w-^-0  agrees  with  the  full  three-dimensional 
theory  by  Lee  and  Chang  (1979). 

Symmetric  plane  strain  wave  propagation  on  the  basis  of  the  theory 

characterized  by  (7.3)  to  (7.7)  yields  equations  for  u2,Uj,uj,Sj,6j,6j.  A 

non-zero  solution  of  these  equations  gives  a  quartic  dispersion  equation  for 
2  2 

5  in  terms  of  o>  .  We  omit  the  details  but  note  that  one  root  of  this  quartic 


is  such  that 
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as  ai-*-0.  This  agrees  with  the  value  predicted  for  symmetric  plane  strain 
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waves  from  the  effective  modulus  theory  of  section  5  which  has  no  dispersion 
relation.  In  the  case  of  isotropic  plates  the  limiting  value  (8.5)  as  oj-*-0 
agrees  with  that  found  from  the  three-dimensional  theory  by  Lee  and  Chang 
(1979) . . 


In  a  similar  way,  the  group  of  equations  (7.8)  to  (7.12)  separate  into 
two  further  groups,  one  for  antisymmetric  SH  waves  and  the  other  for  anti¬ 
symmetric  plane  strain  waves.  In  the  first  group  we  have  nonzero  values  for 

^  2  ’  ^  2  ’  U2  t*ie  seconc*  group  SpUj.S^.Sj.u^  ,u^  are  nonzero.  For  anti- 

2 

symmetric  SH  waves  we  find  a  quadratic  dispersion  relation  for  £  in  terms  of 
2 

to  .  There  is  no  wave  of  this  type  to  be  found  from  the  effective  modulus 

theory  of  section  5.  Also,  in  the  case  of  antisymmetric  plane  strain  waves 

2  2 

there  is  a  quartic  dispersion  equation  for  £  in  terms  of  0)  .  One  root  of 
this  quartic  is  such  that 


st 

2 

co 


p2h2  *  2Plhj 


(V1)2 

1  ,  3r2  11  1  C33J 


,1  11.2 

2.  rl  11  (  c33) 


(8.6) 


X  ,  Jrt  il  JJ  1  ,2  ,3  ,2,  ,,  ,  ,4,  rl  XX  JJ  ,  1 

12  V  C11  ‘  2  33  }  (3  hl  hlh2  +  *  hlh2){  C11  *  1  33 


33  33 

as  This  agrees  with  the  value  predicted  for  antisymmetric  plane  strain 

waves  from  the  effective  modulus  theory  of  section  5  which  has  no  dispersion 

relation.  In  the  case  of  isotropic  plates  the  limiting  value  (8.6)  as  o>-*-  0 

agrees  with  that  found  from  the  three-dimensional  theory  by  Lee  and  Chang 


(1979) . 
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9.  Large  number  of  laminates 

The  theory  of  section  4  and  section  7  may  also  be  used  to  discuss 
periodic  stress  distribution  in  a  stratified  medium  consisting  of  a  large 
number  of  alternating,  parallel  layers  of  two  homogeneous  orthotropic 
materials  whose  directions  of  orthotropy  coincide.  We  consider  only  the 
case  in  which  deformation  and  temperature  fields  are  identical  in  every 
other  layer  so  that 


h.  =  h  ,u=u  ,  6  -  =  6  ,  0  o  =  0  ,<j>=4> 

n+2  n  ~n+2  ~n  ~n+2  ~n  n+2  n  n+2  n 


(9.1) 


In  view  of  the  continuity  conditions  (4.3),  it  follows  that 


u  .  =  u  ,  h  ,5  ,  =  -h  5  , 

~n+l  ~n  n+l~n+l  n~n 

9  ,  =  0  ,  h  ,(J>  -  =  —  h  <+> 

n+1  n  n+1  n+1  n  n 

We  may  now  select  three  typical  laminates,  say  1,2,3,  and  use  the  theory  of 
section  7  in  which 


(9.2) 


u. 

l 


5.  * 

l 


9*0 


=  0 


+ 

u. 

l 


2“i  •  hlSl  ■  -  2h2Si 


0+  =  202  ,  h1$+  =  -  2h24>2 


t 

u~ 


3 


k 

u 


*3 


(9.3) 


The  two  groups  of  equations  (7.3)  to  (7.7)  and  (7.8)  to  (7.12)  simplify 
considerably.  We  restrict  attention  here  only  to  the  group  (7.3)  to  (7.7) 
which  determine  extensional  deformations  parallel  to  the  laminates.  Under 
isothermal  conditions,  (7.3)  and  (7.5)  reduce  to 


V 
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i,2  .2  1  -2 

N  o  +  t.  -  tQ  =  p_u 
a3,a  u  3  u  3  2 


<3,3  '  ^  *  Vu*3  V3>  ■  V2«3 
Na&,a  * 2 -u^e  ‘  uV  *  2piue  - 


(9.4) 


<3,3  -  V*3  *  hl  („l3  *  »*!>  ■  -  iCPAW*! 


and  this  leads  to 


<2hi<3  ‘  h2<3>  ,3  ‘  ‘’Vi  ‘  V?  '  2lh1p2“2*h2pl“l/hl'i3 


‘2N3B*N36),3- 

•  V?  -h. 


(9.5) 


Also,  from  (7.4)  and  (7.6)  we  obtain 


2N2  +  N+ 


,2.11  .  1.1U  2 


,2.11  1.11.  2 


11  T"U  ■  2("A11  +  An^en  +  2(  A22  +  A22)e22  +  2(2ai1  '  h?  All)03 


9m2  +  _  ,,2.11  1.11.  2  , ,2.22  1.22.  2  _r2.  h2  1A  ..2 

2N22  N22  "  2(  A22  A22)ell  2("  A22  A225  22  2 C  A22  '  hj  A22)63 

2  +  2  +  2  12  1  12  2 
2N12  +  N12  =  2N21  +  N21  "  4(  A12+  AiPe12  ’ 

h2 

2hlM13-h2MI3  =  2^hl  2Ell+hf  1e11)P31  ’ 

h2 

2hlM23 '  h2M23  =  2(hl  E22  +  h^  E22^P32  ’ 

2hlV3  ‘  h2V3  =  2^hl  A11  "  h2  All)ell 


(9.6) 


+  2<hl  2A22-h2 


It  is  now  a  simple  matter  to  discuss  harmonic  wave  propagation  in  the 

2 

x- direct ion.  There  is  a  symmetrical  SH  wave  for  which  is  non- zero  and 
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(9.7) 


pl  +  p2 

2A12  ♦  1A12 
*12  ♦  *12 


Plhl  *  P2h2 
,  2  12  ,  1  12 
h2  c12  +  hl  C12 


When  the  layers  are  isotropic,  this  reduces  to  the  result  given  by  Sun, 

Achenbach  and  Herrmann  (1968).  In  addition,  equations  (9.4)  to  (9.6)  yield 

2  2  • 

longitudinal  waves  for  which  u^.S^ 

2  2  2 
2  ^21  2  ^2Piai  22  ^2  1 

I  All  +  Hr  EU)5  -  fp2°2  +  *  A  +  i  A1 

hl  hl  hl 

„  r  .2. 11  1. ll._2  t  .  2.  .2  h2  1.  .2-2  _ 

K  AH+  AuH  '  (p!  +  P2)w  1  -  C  An  -  jp  au)  5=0, 


or 


,1  ,,  2  13  .  1  13..  2  -2  1  *  *.  ..2  2 

*■12  (h2  c13  hl  c13  h2  ^  "  12  (P2h2  Plhl^h2 


n2  ,,  2  33  .  1  33..  v  r,,  2  11  .  1  11.-2  .  *  *  .2. 

+  h^(hl  C33  +  h2  C33^  *  *(h2  Cll+hl  C11H  -  (P2h2  +  P1h1)o)  ] 

f2  11  1  11.2.2-2  _  n 
-  (  c33 "  C333  h2^  -  0  ‘ 


(9.8) 


For  the  stratified  medium  considered  here,  we  use  the  value  8  =  1/12  in  the 

-  2  i 

expressions  for  E^,  in  (2.19).  Again,  when  the  layers  are 

isotropic,  (9.8)  reduces  to  the  result  given  by  Sun,  Achenbach  and 
Herrmann  (1968) . 
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Part  B.  Laminated  Layered  Composite  Rods 


10.  Basic  equations  of  the  linear  theory  of  orthotropic  elastic  rods. 


Consider  a  straight  rod  of  constant  density  p  and  .at  constant 


temperature  9,  bounded  by  the  surface 


F(x,y)  =  0  , 


where  =  (x,y,z)  are  rectangular  Cartesian  coordinates.  The  rod  has 
orthotropic  symmetry  with  respect  to  the  x^-axes  and  has  geometric  symmetry 
with  respect  to  the  x  and  y-axes.  The  basic  equations  of  the  linear  theory 
for  thermoelastic  deformation  of  such  a  rod  have  been  given  in  a  number  of 
previous  papers.  In  particular,  we  refer  to  the  recent  paper  by  Green  and 
Naghdi  (1979b),  which  contains  explicit  values  for  the  constitutive 
coefficients.  The  main  equations  and  constitutive  coefficients  for  iso¬ 
thermal  deformation  of  an  elastic  orthotropic  rod  are  summarized  below  in 
Cartesian  tensor  notation  with  Greek  indices  taking  the  values  1,2  and  Latin 
indicates  the  values  1,2,3. 

We  recall  that  a  Cosserat  curve  51  comprises  of  a  material  curve  £ 

and  a  set  of  two  directors  assigned  to  every  point  of  the  material  curve; 

the  two  directors  are  regarded  as  modelling  the  shape  and  the  deformation  of 

the  cross-section  of  the  rod.  In  the  reference  configuration  of  a  straight 

rod,  let  the  material  curve  be  a  straight  line  £R  which  we  identify  with  the 

line  of  centroid  x  =  y=0  of  the  rod.  Let  u  and  6^  denote,  respectively,  the 

infinitesimal  displacement  vector  and  the  infinitesimal  director  displacement 

of  £r  (or  the  line  x = y = 0) .  Then,  referred  to  the  unit  vectors  e^  along  the 

* 

x^  coordinate  axes,  u  and  6  may  be  expressed  as 
*  — 

The  overbars  on  5  and  6  .  in  Green  and  Naghdi  (1979b)  are  omitted  here  for 
convenience.  ''a  ai 
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u  =  u.e. 


6  =  6  . e. 

~a  ai~i 


and  the  relevant  strain  measures  are 


(10.2) 


Y0S  ■  SaB  *  V  '  ya3  *  Sc.3*3Va2  •  Y33  *  28u3/3z  * 


k  0  =  36  0/3z  ,  <  _  =  3<S  ,/3z 

af3  a&  a3  ct3 


(10.3) 


The  contact  force  vector  n,  the  contact  director  forces  p  ,  the  intrinsic 

director  forces  it  and  the  external  fields  f ,SL  ,  referred  to  the  basis  e., 
~a  -  ~a  ~i 

can  be  expressed  as 


n  =  n.e.  ,  p„  =  p.e.  ,  tt  =  ir.e.  , 

~  i~i  Ect  *011-1  -a.  ai~i 


f  =  f.e.  ,  l  =  l  . e . 
i-i  ~a  ai~i 


(10.4) 


The  basic  equations  of  the  linear  theory  under  discussion  separate  into 
four  groups,  two  for  flexure,  one  for  extension  and  one  for  torsion.  These 
are  summarized  below. 


Flexure  FI 

*23  =  n2  =  k5Y23  *  ml  =  P23  =  k15lc23  ’ 

2  2  2 

3n2/3z  +  pf2  =  p3“u2/3t  ,  -  n2  +  p£23  =  pa23  S^/St 

Flexure  F2 

*13  =  nl  “  k6Y13  ’  m2  =  "  P13  =  "  k16<13  * 

3n^/3z  +  pf^  =  p32u1/3t^  ,  3m2/3z  +  n^  -  pi,^  a  -  paj326jj/3t2 


(10.5) 


(10.6) 
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Extension  E 


*11  =  2klYll + k7Y22 + k8Y33  ’  *22  =  k7Yll  + 2k2Y22 + k9Y33  ’ 


n3  =  k8Yll  +  k9Y22  +  2k3Y33 


P11  =  k10ICll  +  ;*  k17K22  ’  P22  =  ^  k17Kll  +  kllK22 


'  (10.7) 


9n3/3z  +  pf3  =  p92u3/9t2  , 

Sp^/az  +  pX,xl  -  =  pa182<S11/3t2  ,  3p22/3z  +  p&22  -  tt22  =  pa232<522/3t2 


Torsion  T 


*12  =  *21  =  *  k4CY12+Y21}  '  P12‘P21=m3 


P12  =  k12K12  +  ^  k14K21  *  P21  =  kl3l<21  +  ^  k14IC12  ’ 


(10.8) 


3p12//^t  ”  *12  +  P^12  ~  <^12//^{:  ’  ^P21^^  "  ^21  +  p^21  ~  ^*2^  ^21^^ 

We  also  record  here  the  value  of  the  constitutive  coefficients  for  a 
rod  of  rectangular  cross-section  bounded  by  the  planes  x  =  ±h  a,  y  ■ ±  %  h. 

In  terms  of  the  elastic  constants  s^  (or  c^)  for  an  orthotropic  material 

§ 

in  the  three-dimensional  theory,  the  relevant  constitutive  coefficients  are 


*  3  *  3  * 

p  =  ahp  ,  pa^  =  ha  p  / 1 2  ,  pa 2  =  ah  p  /12  , 


(10.9) 


kj=  ^ahcjj  ,  k2  =  J*  ah  c22  ,  k3  =  h  ah  c33  , 


k?  =  h  ah  c22  ,  kg  =  ^  ah  c33  ,  kg  =  \  ah  c22  , 


(10.10) 


,  1  .  3  13  .  1  .  3  23  v  n 

k10  *  12  ha  C13  ’  kll  =  12  811  C23  ’  k17  =  0  ’ 


These  coefficients  are  special  cases  of  those  given  by  Green  and  Naghdi  (1979b) 
but  revised  values  are  given  here  for  k12,k14  and  k4.  We  also  note  that  in  the 
paper  of  Green,  Naghdi  and  Wenner  (1974,  p.  501),  Eq.  (8.55)  should  read 

k12  -  ^ki4  *  /4U  + -jOR  )  ]•  In  addition,  we  now  set  k^sO. 
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0 


(10.11) 


'  k12  =  k13  =  H3)  ’  k4  =  ^ahci2/8 


k15  -  ah5/(12  s”)  ,  k16  -  ha3/(12 


(10.12) 


Previously  y  was  taken  to  be  8.  For  vibration  problems  we  set  y=  5.135, 
the  first  zero  of  the  Bessel  function  .^(S).  The  coefficient  3)  is  the 
classical  torsional  rigidity.  For  a  rod  in  which  a/h  >  (c^/c.^)  2  this 
is  given  approximately  by 

*  •  I  ah3  c»Cl  -  IfL  (!||)\anh[g  (^f|)%]>  . 

ir  a  c*_  c,  _ 


(10.13) 


The  coefficients  k,.,^  are  given  by  expressions  of  the  form  (10.26)  in  Green 

★ 

and  Naghdi  (1979b) .  For  a  rectangular  section  these  are  calculated  to  be 
ha3  _  ,  2  13  lL2  .22  ,1  .  "  (-1)" 


12kA  *  S13  S33  42  2  2_KrMa  7  13.  23,% 


n-1  ir  n  cosh[-~ ^  (sjj/s“)  *] 


22  22  13 

_  2  .13*33  0  4ha  S33  rS13^  ” 

=  ^a  S13(_23  +  2) - r-  (“23}  1 


"i3^  _ m r _ 

23"'  L  “  ,,2  .  ,  r (2n+l)irh  '  23,  13^-, 
s23  n=°  t2n+1)  sinh[-^ — fa —  S23^ S13^  ^ 


(10.14) 


ah3  .  ,/L2  23  2  11  /  1  "  (-l)n 

12k,.  "  2  n  S23  a  S33  42  2  2,  rn7rh  423.13,%,' 

5  n=l  tt  n  cosh[— —  Cs23'si3J  i 

11  ,u  11  23  „ 

_  t,  .2  23  *33  .  4ha  S33  *23,*S  v 


_  j,  2  23  "33-, 

4  h  S23  ~T3  +  2  ' 
S13 


_ r  2  v  _ _ _ 

_2  13;  77 772  .  ‘  r(2n+l)na  ‘  13,  23^, 

tt  s13  n=o  (2n+l)  smh[^  1 


Let  the  external  force  vectors  applied  to  the  plane  surfaces  x  =  k  a,  x  =  -  ^  a 

and  y  =  ^h,  y  =  -%h  of  the  rod  be  denoted,  respectively,  by  t ,  t  and  by 

t ,  „t  and  write 
u— 


The  expression  for  the  coefficients  k£  and  k6  are  recorded  in  two  different 
forms  for  convenience  of  computation  depending  on  the  magnitude  of  the 
quantity  (a/h)  (sj^/s^)^. 
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t  =  t.e.  ,  t  =  t.e.  , 

+~  +  i~i  —  -  i~i 


t  =  t.e.  ,  .t  =  .t.e. 
u~  u  i~i  SL~  Si  x~i 


CIO. 15) 


Then,  in  the  absence  of  the  effect  of  body  forces  in  the  rod,  we  have 
Pf  = 


rk  a  (kb 


p£l 


,  ^  +  J^)dx  + 

-  k  a 

r  k  a 

,  x(uH  +  ji3)dx  + 

-  k  a 


(  t+  t)dy  , 
kb  ~~ 


kb 


^a(+t-_t)dy  , 


(10.16) 


(k  a  r  k h 


p~2  =  ,  ^h(uH '  jt£)dx+  y(+H  +  -^d>r  • 

*  “  ^  3.  '  -  %  h 
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We  suppose  that  a  rod  of  rectangular  cross-section  contains  N  rod-like 
laminated  layers,  each  consisting  of  a  rod  of  rectangular  cross-section  and  of 
homogeneous  orthotropic  elastic  material  at  the  same  constant  temperature. 

A  typical  cylindrical  layer,  the  nth,  has  constant  density  p*  and  is  bounded 
by  the  plane  x  -  ±  -3  a,  y-H^i  ^h^,  2=  ±  %£,  where  a,Hn,hn>£.  are  constants. 
Each  cylindrical  layer  is  fixed  to  its  neighbors  at  its  major  surfaces 
y  =  Hn  ±  ^  hn  so  that 


n+1  n 


?(hn+l+V  ’  n  =  1,2,...,N-1  . 


(11.1) 


The  mechanical  behavior  of  each  cylindrical  layer  is  assumed  to  be 
governed  by  differential  equations  for  rods  in  (10.5)  to  (10.8)  but  with  an 
additional  index  n  attached  to  each  symbol.  Thus,  for  the  nth  layer,  we 
write  n^m^. ...  .k^, ...  ,k^,p  ,p  n>ai>a2>,,'>uin>ut"<  etc.  Since  displace¬ 
ments  and  force  vectors  are  continuous  at  each  interface  it  follows  that 


j.n+1  .n 

$1  =  $1 


n+l  .  ,  ~n+l 
2  '  *  hn.l®2 


u  4h  i 


n~2 


ut11*^1  .  0  (n  =  1,2,.. .  ,N-1) 


(11.2) 


Forces  applied  respectively  to  the  outer  surfaces  y  ■  H,  -  %  h,  and  y  =  H  +%h 

1  1  '  n  n 

of  the  composite  are 


.t1  ,  tN 

U  u~ 


(11.3) 


while  forces  applied  at  the  surfaces  x = h  a  and  x = a  are ,  respectively, 
given  by 


and  hn  <  y  <  Hn  +  Jj  hn  (n  -  1,2,. . . ,N)  .  (11.4) 

The  relevant  equations  of  motion  and  constitutive  equations  for  each 
cylindrical  layer  are  of  the  forms  F1.F2.E  and  T  in  (10.5)  to  (10.9)  with 
an  index  n  added  to  each  variable. 
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We  consider  first  simple  homogeneous  extension  of  the  rod  along  its  axis 


with  zero  applied  forces  on  its  major  surfaces  x=  ±k  a,  y  =  -  h  h^, 

y  =  H.,  +  \  h.,  so  that 
N  N 

+tn  =  0  ,  _tn  =  0  (n  =  1,2, ... ,N)  , 


1 


=  0 


Recalling  (10.2)  and  the  conditions  (11.2)  we  choose  the  components  of  the 
displacement  and  director  displacement  vectors  in  the  forms 


n  _  « n  •  <  pH  |,n 

iij  ~  Lz  f  ^11  ~  ^  9  *22  ^  * 

u2*1  -  “5  ■  *  V/'*1  *  **  . N-‘)  • 


where  l.M.K^u™  are  constants.  The  remaining  components  of  displacements  and 
director  displacements  are  zero.  It  follows  from  equations  (10.7)  for  the 
extensional  case  that  are  constants  and  we  choose  the  constants 

K11  such  that  =  0  for  n*l,2,...,N.  The  equations  of  motion  in  the  form 
recorded  in  (10.7)  are  then  satisfied  if 


where 


n  n*  n  «  n.  n  n 

r  ^  -  0  ,  y  9 


11  "11 


Pn£- 


r  k  a 


/  .  n  .  Hn  j 

(  t_  +  „t,)dx  , 
j-k a  u  0  4  A 


n.n 
P  ln 


\  a 


3L 


\  a 


pn«,"  =  %  h  |  (  t"  -  .t")dx  . 

22  n  J.  %a  u  2  4  2 


In  view  of  (11.2) 3,  (12 . 1) 3  4>  (12.3) x  ,  and  (12.4)^.,  we  have 


(12.1) 


(12.2) 


(12.3) 


(12.4) 


! 


i 


k  a 


-  h  a 


Ut  3  ^ 


0  ,  \ 


a 


J  -  U 


%  a 


t1^  dx  =  0  for  n  =  1 ,2, . , 


..N 


(12.5) 


Since  7t“2  =  0,  it  follows  that 


2k^  Kn  =  -  k"  M  -  kg  L 


(12.6) 


and  hence 


(k?) 2  n  k"k? 

*U  =  [**?  ~  — “ — ]M  +  [2kp  -  — — ] L  , 


1  k” 
2 

k7kQ 

,.n  J_  9 
k 


(ko)2 


n  r o i  n  7  n  ^  it 

n3  =  [2k8  -  — ]M  +  [4k3  "  “h5 rlL  • 

k2 


(12.7) 


With  the  help  of  (11. 2)^  and  (12 . 1) 3  4»  it  follows  from  (12 . 4) 2  and  (12.3)3 
that 


N  N 

2jPn*'l1  =  °  OT  "  °  • 

n*l  n- 1 


(12.8) 


This  equation  determines  M  in  terms  of  L.  Hence,  from  (12.7)  we  obtain 


n  N 
3  _  n 

L  "  2  n3 
n=l 


N 


■  E  [4k, 

n=l 


<^)2 

k!J 


■1  - 


N  k"k"  _ 

{ 21(2ks  -  “=r)} 

n=l  k2 

N  ~  (k?)2 

L  {4k?  -  } 

n-1  1  k? 


(12.9) 


Since  n^  is  the  resultant  force  along  the  rod  and  L  is  the  extension,  this 
gives  the  effective  modulus  for  extension  in  terms  of  the  moduli  of  the 
individual  layers  of  the  composite. 

Next,  we  consider  flexure  of  a  composite  rod  about  the  x-direction  by 
couples  over  the  ends  of  the  rod.  The  solution  of  this  problem  requires 
equations  of  the  form  FI  in  (10.5),  but  because  of  the  conditions  (11.2) 
at  the  interfaces  between  the  layers  of  the  composite,  the  extensional 
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equations  E  must  also  be  used.  We  again  assume  that  the  applied  forces  on 
the  major  surfaces  of  the  composite  rod  are  zero  so  that  the  conditions 

(12.1)  hold  for  the  present  problem.  In  addition,  making  use  of  (10.2)  and 

(11.2) ,  we  assume  the  following  values  for  the  nonzero  components  of  dis¬ 
placement  and  director  displacement  vectors: 

u"  =  -  ^  Pz2  +  Pz  +  Sn  ,  6^2  =  Pz  -  P"  ,  fij  =  M  , 

u*  =  Lz  ,  u^+1  -U3  =  (Hn+1  -Hn)(Pz  -  P)  (n  =  1,2,. . .  ,N-1)  ,  (12.10) 

<S"2  =  Kn  ,  Sn+1  -  Sn  =  %  hn+1Kn+1  +  %  hnK"  (n  =  1,2, .. .  ,N-1)  , 

where  P,P,Sn,L,M,Kn  are  constants.  It  follows  from  equations  of  ‘he  form 
(10.5)  for  flexural  case  FI  that 


n  1  n  n 
mi  =  kisp 


n  n 
n2  ■ 0 


^23  '  0 


f" 

r2 


(12.11) 


Also,  as  in  the  extension  problem  just  discussed,  we  choose  the  constants 
=  0  for  n  =  1,2,. . . ,  N .  Then 


Kn  such  that  rr^ 


n  r  arTl 

"11  =  [4kl  • 


(kn) 2 

=r-i"*  [2ks  • 

k_  k. 


knkn 


knkn  (kn) 

"3  ■  K  -  -sri"  ♦  K  - 

k2  2 

and  equations  (12.3)  must  again  be  satisfied.  It  follows  that  (12.8)2 
still  holds;  and,  in  addition,  we  require  that  there  is  no  resultant  force 
on  the  end  sections  of  the  rod  so  that 


(12.12) 


N  N 

Ei"  =0  ,  in  =0  .  (12.13) 

n=l  11  n=l 


These  equations  express  L  and  M  in  terms  of  P.  As  before,  the  conditions 

(12.3).  -  yield  the  equations  (12.5)  since  zero  tractions  are  specified 
1  ,  / 
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V 


on  the  surfaces  x  =  ±'ja,  h^  h^.  Moreover,  because  of 

(12.5)  the  conditions  (12.11)3  4  are  satisfied. 

The  resultant  couple  on  any  section  of  the  composite  rod  about  the 
x-axis  is 

N 

=  Z  ®"  +  H  nj)  =  DP  ,  (12.14) 

n=l  1 

where  P  is  the  curvature  of  the  line  x  =  y  =  0  introduced  in  (12.10)  and  D  is 
the  flexural  rigidity  of  the  composite  rod.  The  latter  may  be  expressed  in 
terms  of  the  flexural  rigidities  of  each  constituent  layer  of  the  composite  rod 
and  also  other  elastic  coefficients.  We  illustrate  these  results  by  considering 
a  3-layer  composite  rod  in  which  the  outer  layers  1  and  3  of  the  rod  are 
identical  and  flexure  is  about  the  middle  line  of  the  layer  2.  Then, 


h3  *  h1  ,  H2  =  0  ,  H3  =  -  H:  =  %(hx  +  h2) 

*1  Z  _1 

P  =  P  *  k1  *  kj  ,  etc. 


(12.15) 


and  with  the  help  of  (1 2 . 1 2) 2  the  moment -curvature  relation  (12.14) 
simplifies  to 


m 


1 


_  1  2  u  ,  3  K 

2m1  +  nij  -  ^  (n3  -  n3) 

{2klS*k15*2Hl(4R3  '  “jT_),P 

k2 


(12.16) 


Flexure  about  the  y-axis  is  simpler  than  that  about  the  x-axis.  When 
flexure  is  about  the  y-axis,  the  only  nonzero  components  of  displacement  and 
director  displacements  are 


2  - 
Qz  +  Qz 


13 


=  QZ-Q 


(12.17) 


so  that  from  appropriate  equations  of  the  forms  (10.6)  we  obtain 


54. 


1 


■  "l"0  '  f?-°  •  l13S 


(12.18) 


In  view  of  the  surface  conditions  on  the  composite  rod  and  the  continuity 
conditions  (11.2),,  the  last  two  equations  in  (12.18)  lead  to 


-  \  a 
^  %  a 


t?  dx  =  0 
u  1 


%  a 


-  h  a 


x  t"  •  0 
u  3 


(12.19) 


N 


The  resultant  couple  about  the  y-axis  is 

m2  s  V^-Q  ^e  • 

n=l  n=l 

Finally  we  deal  with  the  torsion  of  the  rod  which  is  free  from  applied 
tractions  over  its  major  surfaces.  In  this  case,  we  have  nonzero  components 
of  displacement  and  director  displacements  of  the  form 


(12.20) 


412  --S21  =  82 


u^  =  Fz  -  (Hn  -HjJgz 


(12.21) 


where  8,N,F  are  constants.  These  satisfy  the  continuity  conditions  (11.2) 
and,  from  equations  of  the  form  recorded  in  (10.6)  and  (10.8),  we  have 


]»"  •  0  ,  n"  >  k"[N.F-  , 


n  n  a  n  n  •  n  ^ 

*12  "  *21  =  °  ’  m3  ^2k12  *  k14  ®  =  ®  S  » 


on  =  on  =o 
12  21  ’ 


(12.22) 


where 
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pnfI  =  ,  Cutni  +  lZl)dx  ’  pni12  "  ,  xCut2  +  *t2)dx 


A11?11  -  >- 

p  Z21  “  '* 


(12.23) 


In  view  of  the  conditions  (11.2)^,  (12. 1)^  4>  (12.22)^  7  g,  from 
(12.23)1  2  3  we  ^ave 


tV  dx  =  0  , 

u  1 


t_  xdx  =  0 
u  2 


(12.24) 


Also  by  (12.1),  4  and  (11. 2)^  we  conclude  from  (12. 23)^  and  (12.22)  that 

N  N 

Z  pnln  =  0  or  In  =0  , 

n=l  n=l  1 


(12.25) 


so  that  there  is  no  resultant  force  over  any  section  of  the  composite  rod. 
The  resultant  couple  acting  on  the  composite  rod  about  the  z-axis  is  then 
given  by 


m3*  Vm3-Hnn? 

n=i 

and  the  constant  N + F  is  determined  in  terms  of  8  by  (12.25)2* 

When  the  composite  rod  consists  of  3  layers  in  which  the  outer  layers 

1  and  3  of  the  rod  are  identical  and  flexure  is  about  the  line  of  symmetry 

of  the  middle  layer,  then  the  conditions  (12.15)  apply  and  the  constant  F 
2 

is  chosen  so  that  ^  =  0.  Using  also  (12.25)2»  we  obtain 


(12.26) 


F  +  Hjg  *  0  ,  N  *  0  ,  n^  =-nJ  =  k^B  ,  n^  =  0 


(12.27) 


The  torsional  rigidity  of  the  composite  rod  can  be  obtained  from  (12.22), 

(12.26)  and  (12.27)  and  is  given  by 

m3/6  -  2  01  ♦  S>2  ♦  2kjHj  .  (12.28) 
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It  is  difficult  to  discuss  more  complex  problems  analytically  without 


specifying  the  number  N  of  the  composite  layers,  but  we  leave  the  dis¬ 
cussion  of  static  problems  at  this  point. 
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13.  A  constrained  theory  of  rods 

We  consider  here  application  of  a  constrained  theory  of  rods  which 
has  a  simpler  structure  in  comparison  with  the  more  general  theory  summarized 
in  section  10.  A  derivation  of  this  constrained  theory  is  contained  in  a 
paper  of  Green  and  Laws  (1973).  For  a  linear  elastic  rod  which  has 
symmetries  of  the  type  discussed  in  section  10,  the  equations  again 
separate  into  four  groups.  The  flexure  groups  are  the  same  as  FI, 2  in 
(10.5)  and  (10.6)  but  the  extension  and  torsion  groups  E  and  T  are 
replaced  by  the  simpler  groups  E  jT^.  Thus 


Extension  E 


n3  =  kY33  '  Y33  =  23u3/3z  * 


3n3/3z  +  pf3  =  p32u3/3t2  . 


(13.1) 


Torsion  T 


612  =  ~  521  ’  m3  5612 


3m3/3z  +  p(&12  -  £21)  =  PO^  +a2)326i2/3t2 


(13.2) 


A  further  constraint  is  often  applied  to  the  flexure  part  of  the 
problem.  Thus,  in  (10.5)  and  (10.6),  -*■  0  while  kg  and  k6-*-«°,  so  that 

the  resultant  force  components  nj,n2  are  indeterminate,  i.e.,  not  given  by 
a  constitutive  equation.  Then,  neglecting  also  the  rotatory  inertia  terms, 
(10.5)  and  (10.6)  are  replaced  by 

Flexure  F  1 


2  2 

m^,  =  -  k133  u^/Sz  , 


2  2 

3n2/3z  +  pf2  a  p3  u2/3t  ,  Sn^/Sz  -  n2  +  pj&23  *  0  , 


(13.3) 
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Flexure  F  2 


m2  =  ^^S^/Sz2 


anj/Sz  +  pfj  =  p32u1/3t2  ,  3m2/3z  +  nj  -  pJ-13  =  0  . 


(13.4) 


The  constitutive  coefficients  ^5 >^6  f°r  a  homogeneous  rod  of  rectangular 
section  are  the  same  as  those  given  in  section  10,  but  the  coefficient  k  has 
the  value 


k  =  ?[k, 


klk9*k2k8-k7k8k9]  .  ah 


4kjk2  -  k7 


2s.^ 

j3 


(13.5) 


Also  Q  is  equal  to  the  torsional  rigidity  3)  given  in  (10.13). 

For  a  composite  rod  of  the  type  discussed  in  section  12,  the  theory  of 
the  present  section  may  be  used  if  we  compute  the  coefficients 
and  k  by  comparing  simple  solutions  with  the  corresponding  four  simple 
solutions  for  extension,  flexure  and  torsion  discussed  in  section  12. 
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14.  Harmonic  wave  propagation  along  the  rod 


The  propagation  of  simple  harmonic  waves  along  a  composite  rod  of  N 
constituent  layers  may  be  discussed  on  the  basis  of  the  theory  of  section 
11.  In  general  the  waves  can  be  categorized  into  two  groups,  one  described 
by  equations  of  the  type  E  and  FI  and  the  other  by  equations  of  the  type  T 
and  F2.  To  make  the  analysis  definite.it  is  necessary  to  specify  the 
number  N  of  the  composite  layers  and  we  do  not  discuss  this  further  here. 
There  is,  however,  one  subgroup  of  waves,  purely  of  the  type  F2,  which  can 
be  discussed  for  the  general  case  of  N  constituent  layers  and  we  consider 
this  briefly.  We  assume  that  all  variables  contain  a  factor  exp  i(£z-oit) 
where  £,co  are  constants,  and  this  factor  is  removed  throughout.  Then, 
recalling  the  continuity  conditions  (11-2),  we  choose 


=  <5  ,  uj  =  u  (n  =  1.....N) 


(14. 1) 


and  take  all  the  remaining  components  of  displacements  and  director  dis¬ 
placements  to  be  zero.  Then,  from  (10.6)  and  section  11  we  have 


n_2  n  2,  .  ,  n  , .  n .ui 

(k6?  -  p  <o  )u  -  i  k&  56  =  p  fj  , 


, ,  n  _2  n  n  2  ,  n. .  .  ,  n  *»  n.n 

(k16S  -  p  ♦  k&)6  +  l  k6  £u  =  p  Z 
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(14.2) 


Assuming  that  the  major  surfaces  of  the  composite  rod  of  rectangular  cross- 
section  are  free  from  applied  forces,  by  (10.16)  the  quantities  f^  and 
are  given  by 


H  a 


)  - 


\  a 


Vi  Vi,d* 


n.n 

p 


•  h  a 

-  H  a 


♦ 


4t3)xdx 


(14.3) 


In  view  of  (11. 2)^,  the  surface  conditions  (12. 1)^  ^  and  (14.3)  we  conclude 
that 
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0 


(14.4) 


N 


N 


z  pnf?  =  o  ,  z  pnan 

_ 1  A  _ m 


.  ,13 

n=l  n=l 

With  the  help  of  (14.4)  we  may  eliminate  u,<$  from  (14.2)  to  obtain  the 

dispersion  relation 


c^p2^4  -  po)2  (k6  +  k16£2  +  ai^6?2)  +  k6k16^  =  0  ’ 


where 


__  N  _  N  _  N  _ N 

k6'£,k«  '  ’  P'„V  •  ' 

n=l  n=l  n=l  n=l 

2  —  4  — 

When  £  -*■  0,  one  value  of  ai  -*-k.,£  /p  and  this  is  the  single  value  which  would 

lo 

be  obtained  by  using  the  F£2  flexure  equations  (13.3)  for  the  same  problem. 
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